Stochastic theory of spin-transfer oscillator linewidths 
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We present a stochastic theory of hnewidths for magnetization oscillations in spin- valve structures driven by 
spin-polarized currents. Starting from a nonlinear oscillator model derived from spin-wave theory, we derive 
Langevin equations for amplitude and phase fluctuations due to the presence of thermal noise. We find that 
the spectral linewidths are inversely proportional to the spin-wave intensities with a lower bound that is de- 
termined purely by modulations in the oscillation frequencies. Reasonable quantitative agreement with recent 
experimental results from spin-valve nanopillars is demonstrated. 
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I. INTRODUCTION 

Since the theoretical predictions of Berger' and Slon- 
czewski2ii concerning the influence of a spin-polarized cur- 
rent on magnetization at a mesoscopic scale, many issues con- 
cerning the interplay between spin transport and spin dynam- 
ics in metaUic heterostructures have been brought to light by 
extensive experimental and theoretical work. The primary ef- 
fect on these length scales is a transfer of spin angular momen- 
tum between the current and the magnetization, which may 
lead to magnetization reversal, the generation of spin waves, 
or both. One can account for the spin-polarized current in the 
equations of motion of magnetization M, to a reasonable ap- 
proximation, by an additional torque of the form 



dM 
~dt 



Ms 



M X (M X p), 



(1) 



where / is the current density, Ms is the saturation mag- 
netization, p is the unit vector along the direction of spin- 
polarization, and (T is a parameter that measures the spin- 
transfer efficiency. For particular orientations of magnetiza- 
tion, current polarization and applied fields, such a term can 
act as a negative Landau-Lifshitz form of magnetic damping. 
As such, the existence of stable precessional states of mag- 
netization can be made possible by tuning the applied cur- 
rent such that it compensates the intrinsic magnetic damp- 
ing (averaged over a precession period). These oscillations 
can be detected experimentally through the oscillating volt- 
ages arising from the time-varying magnetoresistance in a 
spin-valve structure. Evidence for such dynamical states have 
been reported for a number of geometries and magnetic ma- 
tei-ials,ViVffl"'i"i'^i'3i'^i'Si'^i'"i'^.'^ which seems to indicate 
that such stable oscillations are a general feature of spin- 
momentum transfer. 

While the frequencies of these oscillations appear to be 
well-understood from spin-wave theorjs^Si^ and numerical 
simulation|22iS2i2^i2L2i there is less consensus concerning 
the origin of the spectral linewidths. Part of the reason for 
this discrepancy may be due to the different experimental 
conditions, materials, and heterostructure geometries used to 
study these oscillations. For instance, it is found that spec- 
tral linewidths for point-contact geometriesSiiiiii are gener- 
ally narrower than those reported for nanopillars,"24i^ although 



recent experimental evidence from more complex nanopillar 
stacks indicates that this is not a general rule.'^'^** Moreover, 
simple macrospin calculations-'* and more elaborate micro- 
magnetics simulations^^'^'' consistently overestimate the spec- 
tral linewidths by at least an order of magnitude. 

In this article, we present a stochastic theory of spin- 
transfer oscillator linewidths based on spin-wave theory. We 
derive a linearized Langevin equation to show that fluctuations 
about the dynamic state, due to thermal noise, gives rise to a 
spectral linewidth that is inversely proportional to the spin- 
wave mode intensity with a lower bound determined by mod- 
ulations in the oscillation frequency. The main features of the 
theory are shown to give reasonable quantitative agreement 
with recent experiments on spin-valve nanopillars. 



II. STOCHASTIC MODEL 

We believe that the current-induced magnetization oscil- 
lations observed in experiment coiTespond to parametrically 
excited spin- wave modes. Parametric excitation here differs 
from more well-known means such as parallel pumping, ^^ ■^ii 
in which parametric resonance is achieved with an oscillat- 
ing magnetic field applied parallel to the magnetization. In 
the case of spin-transfer, the torque in Q allows a means of 
modifying the spin-wave damping, so parametric excitation 
occurs as the resultant damping for a given spin-wave mode 
becomes negative. The exact profiles of the excited modes 
are geometry-dependent, but nevertheless one can derive from 
spin-wave theory^"'^''^' a generic oscillator equation of the 
form (here we follow the notation of Ref. 21 ), 

bk + {iiUko + rk)bk=Tki{b},H,I,...) + fk{t), (2) 

where bk{t) is a complex dimensionless variable (|6| < 1) 
for a spin-wave mode k, Uko the spin-wave frequency, and 
Tk is a phenomenological relaxation rate that encompasses all 
possible damping mechanisms for the spin wave. All infor- 
mation concerning the relevant magnetic energies (Zeeman, 
exchange, magnetocrystalline anisotropy, dipole-dipole, etc.) 
are contained in ujko through the dispersion relation for the 
spin waves. The bk variables diagonalize the linear part of 
the original magnetic Hamiltonian Ho (without current terms) 
and can be found from the usual series of Holstein-Primakoff 
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and Bogoliubov transformations (a discussion pertinent to thin 
films can be found in Ref..,3Z). For finite magnetic elements of 
arbitrary shape, it is possible to find the corresponding eigen- 
modes numerically^'^-^'*-^^ or through the use of approximate 
spin-wave theories. •^^■•^^••^'^ The quantity Tk on the right-hand 
side represents nonlinear terms that include higher-order inter- 
actions and sources of parametric excitations, such as Eq. [2 
We assume that the experimentally observed oscillation cor- 
responds to a parametrically excited spin-wave mode k that is 
largely populated in comparison to the thermal occupation of 
other modes. This mode is taken to be immersed in a thermal 
bath whose effect is represented by the additional stochastic 
force fk{t) acting on bk, where fk{t) is taken to represent a 
white noise associated with thermal magnons, 

(/,WAMt')> = 2rfc^^4,fc'<5(t-0- (3) 

no(ijjfc) is the thermal occupation of magnons and is normal- 
ized by NS — MsV / {giis), the total spin in the system, in 
the same manner as the variables bk- Here, Mg is the satura- 
tion magnetization, V the volume of the magnetic element, g 
the gyromagnetic factor, and the Bohr magneton. 

In the absence of the nonlinear term, it is straightforward 
to obtain from direct integration of Eq. |2l 

Thus, fluctuations in the magnon population rife — {b*fjDk) re- 
lax towards = no{uJk)/ {NS) with a characteristic time 
r^^. At thermal equilibrium, in the absence of driving terms, 
the occupation of thermal magnons must go over to the usual 
Bose-Einstein distribution, uq = {exp[hujk/kBT] — 1)^^. 
For spin-waves in the GHz range, it suffices to use uq = 
ksT/ (hulk)- In what follows, we shall assume that this rela- 
tionship holds for all levels of driving and nonlinearity consid- 
ered here. We therefore neglect the role of nonlinear damping 
on the thermal noise of the system.^^; appropriate estimates of 
such higher-order effects, obtained for parallel pumping, are 
given in Ref.l40l 

As discussed above, a sufficiently large spin-polarized cur- 
rent can result in a stable dynamic oscillation of magnetiza- 
tion. By expanding Eq. ^in the appropriate spin- wave vari- 
ables, neglecting all three-wave interaction terms and retain- 
ing only four-wave terms that lead to self-interaction of the 
single excited spin-wave mode, one finds an oscillator equa- 
tion with the additional stochastic term of the form^ 

bk + (iujk + rk)bk = - \bk\^)bk + h{t), (5) 

where ujk = Lom + Nk\bk\'^ contains a nonlinear frequency 
shift that can become important for large spin-wave ampli- 
tudes. In the absence of the stochastic term, Eq.|5] describes 
the dynamics of a nonlinear oscillator that is capable of sup- 
porting self-sustained oscillations with a frequency ujk- In- 
deed, at the threshold al = Tk with fk = 0, a Hopf bi- 
furcation takes place and a stable orbit with |6fcp — d = 
[a I — Tk)/(jl is made possible for a I > Tk- 

It is interesting to note that the presence of thermal noise 
suppresses completely the current threshold for self-sustained 
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FIG. 1: Thermal dependence of parametric magnon population Up, 
as a function of supercriticality The shaded region corresponds to 
the zero temperature case Upo. The difference between finite- and 
zero-temperature populations is shown in the inset. 

oscillations. To see this, we introduce the nonlinear terms in 
(|5} into © and take the steady-state limit to obtain 

Cnl + (1 - C)nk - n'a = 0, (6) 

where C = al /Tk is the zero-temperature supercriticality pa- 
rameter. It is straightforward to show that the presence of a 
non-zero nj, term leads to a finite parametric magnon popula- 
tion Up = Uk ~ ti'q for all currents / > 0, as shown in Fig. [2 
For sufficiently large currents such that C, is far from 1, we see 
that the is very little difference between the zero-temperature 
and finite-temperature parametric magnon populations. As 
such, we will assume that \bk\'^ = d for self-sustained os- 
cillations holds at all temperatures to simplify the proceeding 
analysis, and use of the actual temperature dependence will be 
made later when discussing linewidths. 

In addition to the stochastic force fk{t), noise sources can 
lead to a direct modulation of oscillation frequency u)k itself. 
Such modulations can originate from external field sources or 
other fluctuating magnetic elements in close proximity. An 
estimate of the frequency modulations"*' can be made by con- 
sidering the action of a random thermal field hr{t) superim- 
posed on the external field i/o, H{t) = Hq + hr{t), where 
hr{t) is assumed to be a white noise with zero mean, 

(hrit)) ^ 0; {hr{t)hr{t')) = 2Dr5{t - t')., (7) 

and Dr is an appropriate diffusion constant. The action of the 
random field on the equations of motion for the b variables are 
obtained by supplementing the original diagonal Hamiltonian 
Ho with the random field Hamiltonian 

Hrit) = -fiJ2hrit)S,, ~ hu;r{t)Y,b*kbk, (8) 

i k 

where ^ = g^QfiB, ujr(t) = jhr{t), and 7 = y/io^s/Zi is 
the gyromagnetic constant. Following the usual procedure for 
generating the equations of motion, bk = i[7io + ^r, bk]/fi, 
it is a simple matter to show that the random field leads to a 
modulation of the spin- wave frequency uj{t) = ujq + oj'{t)- 

In experiment, the current-driven oscillations in spin-valve 
devices are detected through a time-varying voltage due to the 
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corresponding magnetoresistance oscillations. Giant magne- 
toresistance is related to the relative orientations of the free m 
and fixed magnetizations AI ~ Mz, R = Ro + AR {M ■ rfi), 
where Rq and Ai? are material-dependent constants that can 
be readily determined in experiment. We allow for a misalign- 
ment between the easy axes of the free and reference layers, 
which is denoted by the angle measured from the z axis in 
the xz plane. The corresponding (scaled) voltage fluctuations 
are then given by 

V = ^ j^^" = m^cosg - ^ sin6'(Tn-|- + to_), (9) 

where m± — nix ± imy, and /o is the average dc current. Ex- 
panding m in the appropriate spin-wave variables and assum- 
ing that only the parametrically excited mode k contributes to 
the time-varying magnetization, we find that the two-time cor- 
relation function for the voltage fluctuations, {v{t)v{{))), can 
be written as 

{v{t)vm ^ F{k,6){{hk{t)him + c.c.) 

+ G{k,e){bu{t)bl{t)bmbm)-. (10) 

where F{k, 6) and G{k, 6) contain all the factors related to 
the transformation to spin-wave variables. To obtain dlOt . 
we have assumed that there exists a center of symmetry such 
that fefe = 5-fe and that the spatial average of hk is non-zero. 
These assumptions are justified in rectangular magnetic ele- 
ments, for example, in which the lowest-lying excitations have 
been shown to be cosinusoidal standing waves in the dipolar- 
dominated regime. In the absence of misalignment {9 — 0), 
contributions to the magnetoresistance oscillations can only 
come from the term, which gives rise to the G(fc, 9) terms 
to lowest order and result in spectral frequencies that are twice 
the oscillation frequencies. In what follows, we assume that 
the misalignment is non-zero and that the two-time correlation 
functions in the spin-wave variables, proportional to F(fc, 9), 
are sufficient to describe the experimental data. 



III. CORRELATION FUNCTIONS 

The Fourier transform of (u(i)w(O)) gives the spectral 
density of voltage oscillations measured in experiment, so 
the main goal of this section is to obtain expressions for 
{bk{t)bl{0)) and (6*(t)6fe(0)) from To begin, we make 
the assumption that fluctuations in the spin-wave amplitudes 
are sufficiently small so as not to modify the frequency in the 
stable oscillation regime. This allows a rotating-wave approx- 
imation [bk{t) — i3fc(t)e^*'^*] to be made to obtain 



Bk^P{dk~\Bk\^)Bk + ^fk{t), 



(11) 



where /3 = al, 4 = {al - Tfe)//?, ^fkit) = e'^'fk{t), 
and q — Tkno/2NS. It is reasonable to assume that the 
stochastic noise remains white in the rotating frame, so we 
will take fk{t) to have zero mean and satisfy {fk{t)fk' (^')) = 
A5k ki5{t — t'). For the sake of clarity and without loss of gen- 
erality, we will omit the mode index k in what follows. Fol- 
lowing the treatment in Ref.|42i we introduce B = bi + 162 



and f = fi + if2 to write Eg. II lias 

k^f3{d^bl~bl)b, + ^Mt), 



(12) 



where i = 1,2 and {fi{t)fj{t)) — 2SijS{t — t'). Because 
we are interested in separating out the amplitude and phase 
fluctuations in B, we find it conv enient to introduce B = 

Aexp{i(j>), where A = \/b\ +63 and = iwc'^ib^lbi), 
which leads to 

i - - = cos + /2 sin0); 



A 



■(-/isin(/) + /2C0S( 



(13) 



Instead of solving ( I13l l. one can equivalently consider 

i-/3(d-A2)A--| = ^/^(t); 



with 



{fA{t)fA{t')) = {Mt)uit')) = m~t'), 

{lA{t)U{t')) = 0, 



(14) 



(15) 



as these Langevin equations lead to the same drift and diffu- 
sion coefficients^^ (in the Stratonovitch definition) and there- 
fore describe the same physical system. 

The stochastic equations MA\ represent a set of coupled 
nonlinear Langevin equations. In general, it is difficult to 
compute correlation functions directly from these equations 
and one can instead solve the corresponding Fokker-Planck 
equation. For our spin-transfer oscillators, we note that for ex- 
citations sufficiently above threshold (determined in this case 
by the thermal noise floor), one can linearize the Langevin 
equations about the stationary state = d and consider 
only the small fluctuations about this dynamic equilibrium. 
Let = \/d + a(t), where a(t)/Vd ^ 1. If we retain only 
linear terms in a{t) and neglect a term proportional Xo q/d, we 
find a simple Ornstein-Uhlenbeck process for the amplitude 
fluctuations and a Wiener process for the phase fluctuations. 



(16) 



where k — 2(3d — 2(<t/ — T). In the absence of the stochastic 
terms, (I16> shows that amplitude fluctuations are damped out 
with a characteristic frequency k that increases as the driving 
current increases (or equivalently, as the system moves fur- 
ther away from threshold). While this can lead to a linewidth 
broadening with the driving current, it is shown later that it is 
the phase fluctuations which dominate in the large amplitude 
regime. 

From the formal solutions of il6\ . it is straightforward to 
evaluate the desired two-time correlation functions for large 
times t,t' ^ after which initial correlations are unimpor- 
tant. We find 

{a{t)a{t')) = ie-''l*-*'l; 

K 

m')) ^ g-ilt-t'l/d^ (17) 
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which shows that phase fluctuations are inversely proportional 
to the mode intensity d and includes a temperature depen- 
dence through q. 

While the last result demonstrates that phase fluctuations in 
B are inversely proportional to the mode intensity d, there 
exists a lower limit to the phase noise arising from modu- 
lations in the oscillation frequency, which are implicitly ex- 
cluded from the analysis above by virtue of the rotating-wave 
approximation. An estimate of this additional contribution to 
the phase noise is found from 



where 



-iu}(t')\ _ gii^oig- 



-7^_Dr|t-t' 



(18) 



which shows that frequency modulations will lead to an 
amplitude-independent linewidth in the spectral density that 
is proportional to j^Dr- If we take a random magnetic field 
corresponding to superparamagnetic fluctuations of the entire 
magnetic element (Goldstone mode),**^ we find 



70; 



(19) 



where a is the Gilbert damping constant. We point out that 
it is also possible to include the effects of spin-shot noise into 
this definition by an appropriate renormalization of Other 
fluctuations in the applied current will not be considered here. 

The spectral density of magnetoresistance oscillations ob- 
served in spin-transfer experiments is proportional to the 
Fourier transform of the correlation functions {b{t)b* (0)) and 
{b* {t)b{0)) . The contribution to positive frequencies is given 
by {b* {t)b{0)) , which in terms of our definitions for A and 
is 

{b*{t)b{0)) = (e*'^We-*'^(o)>(B*(t)B(0)), 

= (e»'^(t)e-''^(o))(^(i)^(0))(e*-^'(*)e-*'^(o)), 

^ K 

where in the second line we have used the fact that amplitude 
and phase fluctuations are uncorrected in the linearized equa- 
tions. Because the correlation function is a sum of two ex- 
ponential functions, the spectral density will be a sum of two 
Lorentzians centered about oj — luo, S{u!) = d Si{lu; Tf^) + 
{q/ k)S2(uj; r|®), where Siioj; T) is a normalized Lorentzian, 
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Tt{uJ- CJo)2 + r2 ' 



(21) 



with Tf = j^Dr + q/d and Vf = Tf + k. While the pres- 
ence of the q/n and k factors from the amplitude correlation 
function appear to introduce a complicated d dependence. 



1 Ad 
K = -1 - 



n'Q{Lo) I ~ d 



2 1-d' K 



(22) 



one can show numerically that, in practice, such terms do not 
significant modify the overall 1 /d dependence of the effective 
linewidth of S{uj), and as such, one can use 



Sold) 



(23) 



5*0 (d) 
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Id' 






r {l-df 
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4^2 - d) 



(24) 



to a good approximation. 



IV. RESULTS AND COMPARISON TO EXPERIMENT 

To obtain the full temperature dependence it is necessary 
to replace d with the correct parametric magnon population at 
finite temperatures. 



np{d,T) = - + .^rf2 + _ d)n'oiT) - n',{T), (25) 

which is obtained from the solution of Thus, the spectral 
linewidth is given by 



np{T)' 



(26) 



This is the main result of this paper All material parameters 
are encompassed in the spin-wave relaxation frequency Tk, 
while the thermal and intensity dependences of the linewidth 
are contained in k-RT/np{T). If we apply the Gilbert model 
for spin-wave relaxation and ignore propagating losses, the 
damping rate can be obtained from.^' 



r 



du! 



(27) 



H 



where ujh represents the circular part of the magnetization 
precession. Thus, the spectral linewidth for a parametrically 
excited spin-wave mode bk with Gilbert damping, subject to 
thermal noise sources Q and ( I19> . results in an appealing 
simple form 



Aujk 



2^a 



1 dujk 1 

2 duH d 



k^T. 



(28) 



This result shows that random field fluctuations given by 
([19} play a minor role compared to the thermal noise, since 
doJk/duJH is of order unity and d <C 1 for most practical sit- 
uations. We will neglect the random field term in what fol- 
lows, but we would like to emphasize that other external field 
sources not considered here could, in principle, give more sig- 
nificant contributions. 

A plot of the linewidth factor k^T/up as a function of d 
is presented in Fig.|2l Here, we have used a parametric spin 
wave frequency of ujko/i^TT) — 12 GHz for the scale factor in 
the temperature dependence. For low excitation amplitudes, 
one observes that the linewidth is nearly independent of the 
intensity d and is instead dominated by the thermal noise term 
jiq. As parametric excitation is increased past d ~ 10^2^ the 
linewidth decreases rapidly like the inverse of the intensity 
1/d. We would like to point out that such a spin-wave theory 
is expected to breakdown somewhere in the region d > 0.2, 
in which it would be necessary to retain higher-order terms in 
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FIG. 2: Linewidth factor keT/up as a function of d for several tem- 
peratures. 
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FIG. 3: Linewidth factor k^T jup as a function of temperature for 
several spin-wave intensities d. 



for a correct description. Furthermore, there is no guaran- 
tee that the single parametric mode will remain stable in this 
highly nonlinear regime. One should therefore regard the pre- 
dictions for d> 0.2 as an illustration of the general qualitative 
behavior. 

The renormalization of the parametric magnon population 
by the thermal magnons in \25\ leads to a non-trivial ther- 
mal dependence for the spectral lines, as shown in Fig. |3l 
In the low-amplitude regime where thermal magnons domi- 
nate the linewidth, one observes that the linewidth increases 
like Vr, due to the corresponding term under the square- 
root in M51 . Such a behavior has already been observed in 
experiment^^ and can be understood from random-field argu- 
ments. As the mode excitation increases, a transition towards 
a linear temperature dependence occurs between d = 0.01 
and d = 0.1. In this regime, the thermal magnons are essen- 
tially separated from the paramagnetic magnons, so their role 
is simply to provide an independent thermal noise term with a 
linear temperature dependence. Such a regime change should 
be detectable in experiment. 

Much of the discussion so far has made use of d as a 
characteristic parameter. Indeed, the results above show that 
quantitative predictions about the spectral linewidths can be 
made without knowledge of the spin-transfer efficiency a (and 
therefore the absolute values of the threshold currents), which 
can only be obtained from detailed transport calculations 4^ 
In the same spirit, we believe that quantitative comparisons to 



Co 



10 



10" 



'10^ 



10 



10^ 





0.001 0.01 0.1 1 1 10 100 

d T(K) 

FIG. 4: Spectral line intensity Soiuj = uoko) as a function of (a) 
intensity d and (b) temperature T. The symbols in (a) are a guide to 
the eye. 



experiment can be made through a characteristic feature of the 
parametric spin wave spectra that is independent of the cuiTent 
and field configurations. We suggest that the intensity of the 
spectral lines, given by MAI , is a good experimental parameter 
against which linewidths can be compared. In Fig.|^, the the- 
oretical spectral line intensity 5*0 in MA\ is shown as a function 
of mode intensity d. As one would expect, the power pumped 
into the parametric mode increases as the mode intensity in- 
creases, with a low-amplitude region that is essentially flat 
due to thermal magnons. Departures from this monotonic in- 
crease, observed in any experiment, would be a sign of mode 
instability and the generation of other spin- wave modes. A 
stronger temperature dependence is seen for large amplitudes 

(Fig.ia^). 

In the regime d > 0.02, the spectral line intensity is approx- 
imately proportional to the square of the spin-wave intensity 
d. Therefore, one can expect the relationship Aujk 1/Sq 
to hold for large excitation amplitudes. An example of this 
behavior has been observed in recent experiments on rectan- 
gular spin- valve nanopillars.^^ The experimental results with 
a comparison of (|2l and ^ for large d are shown Fig. [S] 
The free-layer in the spin-valve stack has dimensions of 100 
X 50 X 2.8 nm, with (IqMs = 0.85 T and a Gilbert damping 
constant of a = 0.02. The data are taken from a low-current 
regime in which only one oscillation mode is observed and 
no evidence of mode instability is detected. The solid lines in 
Fig. |5] are based on the theoretical prediction of ( I28l l. where 
the smallest linewidth value from the 150 K data was used to 
obtain a correspondence between the measured spectral line 
intensity 5*0 and the parametric magnon population Up; this is 
the only fitting parameter. To calculate the spin-wave damp- 
ing rate (Eq. \27l . we used the dispersion relation appropri- 
ate for rectangular magnetic elements and assumed that the 
lowest-order spin-wave mode is excited. All other parame- 
ters are obtained from experiment. One observes that an ex- 
cellent fit is obtained for the lowest temperature data set (150 
K), with progressively larger deviations as the temperature is 
increased, the origins of which are still under investigation. 
Nevertheless, the theoretical curves give a good quantitative 
bound on the experimental linewidths. We note that the ran- 
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FIG. 5: (Color online) C omp arison of Eg. 1281 (solid lines) with ex- 
perimental data from Ref. 28. The only fitting parameter is the scale 
factor relating the measured mode intensity So tx rip to the para- 
metric magnon population rip, which is obtained from one data point 
from the 150 K series. 



dom field contribution il9i would lead to a limiting linewidth 
of about 0. 1 MHz with these experimental parameters. 



V. DISCUSSION AND CONCLUDING REMARKS 

The quantitative agreement between our stochastic model 
and experiment lends credence to the hypothesis of paramet- 
ric spin-wave excitation in spin-transfer oscillators.^"-^' This 
description is most appropriate for nanopillar geometries in 
which the eigenmode spectrum is well defined. For point- 
contact geometries, self-localization of spin-wave excitations 
has been shown to be crucial to describe experimental data."*^ 
One could, in principle, extend the ideas in this article to 
point-contact geometries by using the well-studied stochas- 
tic Ginzburg-Landau equation."*^ The relevance of Brownian- 
motion for the parametric spin waves also suggests why it is 
inherently difficult to reproduce narrow linewidths by straight- 
forward numerical time-integration of magnetization dynam- 
ics. In all time-integration schemes, to the best of our knowl- 
edge, the effect of finite temperatures is accounted for by 
including random noise terms to the equations of motion of 
magnetization, as we have done. While one could expect am- 
plitude fluctuations to reach dynamic equilibrium in a time 
of (e.g., for C = 1.25, ~ 5 ns for the nanopillars 
in Ref. 28), the ensuing phase fluctuations would require in- 
tegration times to at least the order of 1/A/ ~ 0.1 /is for 
meaningful comparisons with experiment. Furthermore, it is 
necessary to average over a large number of realizations of the 
stochastic force to recover the appropriate statistics. This rep- 
resents a computationally-intensive task, so any lighthearted 



venture into such a program of study will necessarily lead to 
inconsistent results. 

Within the framework of the stochastic model presented 
here, spectral linewidths of current-driven magnetization os- 
cillations are due to phase noise that is determined purely 
by the spin-wave intensity and temperature, where the latter 
controls the level of thermal noise due originating from ther- 
mal magnons present in the system. This relationship holds 
independently of any details concerning the spin-transfer ef- 
ficiency of the heterostructure a, notably its magnitude and 
angular dependence. As such, no attempt is made to discuss 
absolute linewidths in terms of any external control parameter 
such as the applied electric current. Furthermore, we expect 
that the excited spin-wave mode will not be stable to all levels 
of current-driving, and therefore details concerning the non- 
linear instabilities would be required to give any quantitative 
predictions on the minimum linewidth attainable for a partic- 
ular system. 

In the discussion of frequency modulations due to external 
sources, we evoked briefly the random fields il9i appropri- 
ate for the superparamagnetic fluctuations of a single-domain 
particle4^ Indeed, transverse components of this random field 
can also give rise to a stochastic term of the form /o(i) in 
where only coupling to the (quasi-)uniform mode (k = 0) is 
possible given that the random field is taken to be uniform 
across the magnetic particle. As for the frequency modula- 
tions, it is easy to show that the transverse components also 
lead to an amplitude-independent phase noise of the order of 
0. 1 MHz. This does not appear to be the dominant mecha- 
nism for spin-valve nanopillars at least, where a strong 1/d- 
dependence is seen. 

In summary, we have developed a stochastic theory of 
spin-oscillator linewidths based on spin-wave theory. We 
have derived linear Langevin equations for fluctuations about 
dynamic equilibrium, which are shown to lead to spectral 
linewidths that are inversely proportional to the spin-wave in- 
tensity. Lower bounds for the linewidth are shown to arise 
from random fields leading to frequency modulations. These 
results are shown to give good quantitative agreement with 
experiments on spin- valve nanopillars. 
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